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We study the Friedmann equation for the warped codimension-two braneworld background which 
most closely resembles the Randall- Sundrum model. Extra matter on the (Planck) 4-brane, with 
equation of state pe = {a — l)p for the azimuthal pressure, is required to satisfy the junction 
conditions. For 1 < a < 5, we show that there are two static solutions to the Einstein equations 
for given values of the brane stress-energies. Close to the static solutions, the relation between 
Hubble expansion rate H and brane tension reproduces the standard 4D result for small H, but 
exhibits unusual deviations when H is of order the AdS curvature scale. The two static branches 
for 1 < Q < 5 are shown to come together smoothly at a maximum value of H; however the radion 
is shown to be unstable in the branch with higher H . This remains true even with a mechanism 
for stabilization of the radion, i.e., the Goldberger-Wise (GW) mechanism, since large enough H 
overcomes the force of stabilization. Even in the unstabilized case, cosmological constraints on the 
time and spatial variation of Newton's constant are typically satisfied; only fifth force constraints 
require the stabilization. For a > 5 the model is intrinsically stable, without the need for a GW field, 
and in this case we show that infiationary predictions can be modified by the nonstandard Friedmann 
equation; in particular it is possible to get an upper limit on the spectral index, large deviations 
from the consistency condition between the tensor spectrum and ratio r, and large running of the 
spectral index even though the slow roll parameters remain small. 

PACS numbers: 



I. INTRODUCTION 



The cosmology of braneworld scenarios has been widely studied, particularly in the simplest case of codimension- 
one branes, involving only a single extra dimension (see [U, [3j S lj| for reviews). It was initially noticed that the 
Friedmann equation was strongly modified from its usual dependence ^ p to the form ~ p , in contradiction 
to big bang nucleosynthesis and other cosmological tests. Before the realization that such pathological behavior was 
linked to the failure to stabilize the extra dimension, it was discovered that in the Randall-Sundrum (RS) model @, 0| 
with a warped extra dimension due to a negative bulk cosmological constant, the Friedmann equation took a more 
interesting form. 
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where T is the brane tension [1, [^, [l^. Thus the exotic linear dependence H ^ p could arise as a high-energy 
correction to the normal low-energy dependence. The behavior ([T]) is only valid for the single-brane version of RS; in 
the case of two branes, it is necessary to stabilize the interbrane separation [llj, which leads to a more complicated 
dependence for the high-energy corrections p^ . Other interesting forms are possible for the Friedmann equation in 
codimension-one brane models, as in the DGP model [T3| . whose Lagrangian includes an extra Einstein-Hilbert term 
localized on the brane, and leads to a quadratic equation determining H, 

= pT l-H, 2 

in terms of the 5D Planck mass Mc, and 4D Planck mass M4. More intricate possibilities also exist, for example 
through the addition of a Gauss-Bonnet term to the DGP model 

There has also been considerable interest in 6D, codimension-two braneworld models, in part motivated by the 
suggestion of a self-tuning mechanism for the cosmological constant which was argued to be a special feature of 
codimension-two branes fl^, [H, [13, [3 , although this claim is controversial [l^, [l^l . Regardless of the cosmological 
constant problem, it is still interesting to consider the predictions of codimension-two branes for the expansion of the 
early universe, since deviations from the general- relativistic prediction can help constrain the models, as well as lead 
to new possibilities for inflation [21], [S^l . 

In the present work we revisit the problem of modifications to the Friedmann equation, focusing on a warped model 
which is the natural extension of the RS model to six spacetime dimensions 23]-[25]. Part of our motivation is a claim 
that the standard Friedmann equation is not recovered even when the extra dimensions are stabilized. We will 



2 



show that this claim was erroneous, and that ahhough the unstabihzed version of the model indeed exhibits unusual 
expansion behavior at high energies, the expected results of general relativity (GR) are recovered at low energy in 
the case where the radion is stabilized through a Goldberger-Wise (GW) [1^ mechanism, and even in the unstabilized 
model. 

The detailed form of the modified Friedmann equation depends on the equation of state a of extra matter which is 
placed on the 4-brane that plays the role of the UV branc in our model. We find that for most values of a, there is 
an exotic branch of the Friedmann equation which has the Hubble rate H increasing as the brane tension decreases. 
Through a fluctuation analysis, we show that the radion is unstable on this branch, even in the version of the model 
which includes the GW mechanism. 

On the conventional branch of the Friedmann relation, we obtain interesting deviations from the predictions of GR 
near (but still below) the point where the radion instability begins. We show that it is possible to obtain potentially 
observable deviations from the consistency condition which relates the tensor spectral index to the tensor-to-scalar 
ratio, in a model of chaotic inflation on the brane. We note that there is a suppression of the spectral index in the 
braneworld inflation model relative to standard chaotic inflation. Furthermore it is possible (though it requires more 
fine tuning) to get large running of the spectral index. 

The organization of this paper is as follows. In sectionlTTl the basic setup is presented. We derive the bulk equation of 
motion and the jump conditions for the codimension two background with the 3-brane undergoing de Sitter expansion. 
We obtain an exact solution in this background. In section IIIIl we study the Friedmann equation in the model in 
the absence of radion stabilization. The Friedmann equation is obtained using an analytical approximation valid at 
low H, as well as numerically for arbitrary H. We find that standard cosmological expansion is recovered at low 
H, but interesting deviations occur when H ^ £, where £ is the curvature scale of the AdS background. In section 
IIVI we repeat the analysis in the more realistic version of the model which includes stabilization of the bulk using 
the Goldberger-Wise mechanism. To check the effectiveness of the stabilization as a function of the Hubble rate, we 
analyze the stability of the model, with and without the GW field, in section|Vl To explore the physical consequences 
of the modified Friedmann equation, we study chaotic inflation on the brane in section IVII Section IVIII gives our 
conclusions, including a summary of the most important results of the paper. In appendix 1X1 details of the derivation 
of the perturbative approximation to the modified Friedmann equation are supplied. Appendix [B] gives an index of 
the many symbols used in this paper as an aid to the bewildered. 



II. THE MODEL 



We will consider a codimension-two braneworld with a negative cosmological constant in the six-dimensional bulk 
spacetime. The geometry describes a warped conical throat, which is bounded at large r, namely r = P, by a 4-brane 
around which orbifold boundary conditions are imposed, analogous to the Planck brane of the RS model. In addition 
there is a 3-brane at the infrared end of the throat, at the position r — g, which can be thought of as the Standard 
Model brane. The action is 
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(V0)2 + Vi<P) 



-.gA-branc- / d'^Xy/^T^, 



(3) 



where A:| is the 6-dimensional gravitational constant and Qab is the 6-dimensional bulk metric. The 6D cosmological 
constant is negative, Ag = — 10/£^, and gives rise to an approximately AdSg bulk geometry with curvature length 
scale i. We denote the induced metrics on the 3-brane and the 4-brane by g^i, and gab, respectively. T3 is the 3-brane 
tension. Here, >C4-bianc is the Lagrangian density of the matter on the 4-brane that includes the 4-brane tension T4, 
as well as some additional component which is necessary for satisfying Israel matching (jump) conditions at r = P. 
The 4-brane at r = P is needed in order to have a compact bulk and localized gravity. 

The 6-dimensional Einstein equation derived by varying the above action with respect to takes the form 



TAB + Ae^AB — feg^As] 



9 Sab S'X S'b ^{r -P)~ klrs^g^, S'^ S^'^ {r - g) 



(4) 



where J^^-* denotes the 2-dimensional delta function with support at the position of the 3-brane, Sab is the 4-brane 
stress-energy tensor (which we will specify below, eq. ()13|) ) and Tab is the bulk stress-energy tensor. 



Tab = 



QABd''<\>doci>-Y{<\>)GAB^ 



(5) 
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The pure gravity model has a modulus, the radion, which must be stabilized to make a realistic model [16|. For 
this reason we have included a bulk scalar field, which can give the radion a mass by the Goldberger-Wise mechanism 
[2^ . provided that (j) couples to the 4-brane in £4-brano- 



A. Bulk equations of motion 



For simplicity we impose azimuthal symmetry on the extra two dimensions and therefore require that the metric 
depends only on the radial coordinate. The line element has the form 

ds^ = a{r) [-df + e^"'6^Jdx'dx^ ] + f{r)K^ dO^ + JJ^^r^^ (6) 

where _ff is a rescaled Hubble parameter,^ and X is a dimensionful parameter which determines the deficit angle at 
the 3-brane, where there is generically a conical singularity. To insure that the position r — g oi the 3-brane represents 
a single point in the extra dimensions, we require that f{g) — 0. The pure gravity model has a modulus, the radion, 
which must be stabilized to make a realistic model. For this reason we have included a bulk scalar field, which can 
give the radion a mass by the Goldberger-Wise mechanism [2^, provided that (/) couples to the 4-branc in £4-branc- 
The Einstein equations for the model are 

^6 , U2 
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(7) 



where / denotes the derivative with respect to r. 

The variation with respect to (j) gives the Klein-Gordon equation. 



□ 0-^ = 0. (8) 



B. Jump conditions 

To completely specify the geometry, we must consider the jump conditions which relate the metric and bulk scalar 
to the sources of stress-energy at the boundaries. First, let us consider the jump condition for the 3-brane. Recalling 
that the position of 3-brane satisfies 

f{Q) = , (9) 

the deficit angle is given by 

(10) 

r—Q 

A nonvanishing deficit angle indicates the presence of a conical singularity at the position of the 3-brane, related to 
the brane tension T3 by A6' — k^Ts. The jump condition a,t r — g can thus be written as 

K^^(l-kl^). (11) 



A0 = 2tt 



K 



fir) 



^ Since a{g) ^ 1 at the 3-brane, t is not the proper time. The Hubble rate with respect to the proper time is H = H / \/a{g), hence our 
distinction between H and H. 
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In what follows, it will not be necessary to include any coupling of the scalar field to the 3-brane, so the boundary 
condition of (/) at r = g is simply 



- 

Next, for the 4-brane, we define the stress-energy tensor of the 4 -brflllG Sab 



(12) 



T4 



- [Ti + {1- a) — \ gee 



(13) 



The parameter a, which allows for a difference between the 4D and the 60 components of the 4-brane stress tensor, must 
be nonzero to accommodate a static solution with P < oo. The physics which could account for the T4 contribution to 
the stress tensor can arise in several ways [17, 27]. Conservation of Sat implies that the energy density of the source 
T4 must scale like l/i" [H, [3] where 2ttL4 is the circumference of the 4-brane. If T4 is due to the Casimir effect of 
massless fields living on the 4-brane, then a = 5, whereas if it is due to smearing a 3-brane around the 4-brane then 
a = 1. 

We assume Z2 orbifold boundary conditions at the 4-brane, so that the radial derivatives of the metric components 
change sign as one crosses the 4-brane, and their discontinuity is determined by the 4-brane stress energy components. 
Thus the jump conditions at r = P are given by 



/ 
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2fK'^] ^kl^SeeS{r~P), 



(14) 



(15) 



where the term ^/f in the r.h.s. of each equation comes from the ratio of the determinant of the 4-brane induced 
metric to the bulk metric, \J —gj \/ —Q ■ 

By integrating the above equations across r = P, assuming Z2 orbifold symmetry, we get 
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To stabilize the model with the bulk scalar field, we will allow for the possibility that 
through a potential Vp(0) contained in £4_brane- This leads to the boundary condition 

^, 1 dVp 



2f{P) 



(16) 

couples to the 4-brane 
(17) 



at r — P, using the assumed Z2 orbifold symmetry. 



C. Background solution 

The general solution to Einstein's equations for the metric ^ is 

The constants of integration vq, ri can be set to convenient values by rescaling coordinates, — > Ax'^, r — > Br. 
Under this tranformation, the metric functions change as 

a-^aA^B^, H^H/A, f-^f/B"^, K ^ KB (19) 

By choosing A = {^Y^^ and B = {r^r^Y/^ , we can thus set ro ~ ri ~ i, for convenience. To simplify the equations 
however, we will express r in units of i henceforth, so that r (as well as g, P) becomes dimensionless. In the limit of 
H ^ 0, the metric reduces to the AdS soliton solution [28j . 
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Since we have chosen to normahze a{g) = at the 3-brane instead of the usual value, a = 1 the observed Hubble 
rate on the 3-brane is given by H = H / g, and the position of the 3-brane, defined by eq. ©, satisfies 

1 + H^e^ = 4 , (20) 
g^ 

The jump condition for the 3-brane, eq. (|f f p . is then 



and the jump conditions for the 4-brane, eqs. fTCjl . can be expressed as 



v/P2 -p-3 + ij2^2^2 a-1 P 2{a~\) ^ 
^ P+lP-* _ ^P2 _ p-3 + g2g2^2 ^ £fciar4 

^p2 _p-3 +^2^2^2 P 2if"(F2_p-3+^2g,2^2)a/2 > 

These conditions cannot be solved analytically for general values of the Hubble rate. However, we can solve them 
when TJ = 0, as we next demonstrate. 



D. Static solutions 

We wish to find the value of the 4-brane position P and the critical 3-brane tension which satisfy the jump 
conditions for a static geometry, with H = 0: 

P+IP-^ , 3a + lVW-p-^ kla ^^^^ 



^p2 „p-3 a-1 P 2{a-l) 



P + |p-4 - P-3 _ iklar^ 



(25) 



^p2 _p-3 P 2ii'°(P2 -P-3)"/2 

At first sight these equations appear intractable for an analytic solution, but if we first square both sides of eq. (I24|) . 
we find that it becomes quadratic in P^: 

(64 - c)a^P^° + (ca^ - 16a(3a + 5)) P^ + (3a + 5)^ = (26) 

where we defined c — (fc|i?r4)2. The solutions are given by 

From the discriminant A of the quadratic equation, it is necessary that c > Cmin = 60-f 40/a— lOO/a^ = 64— (fO/a— 2)^ 
to have real solutions for P. Furthermore, physical solutions must have P > g = f, where from eq. (|20p the static 
position of the 3-brane is ^ = f . We find that it is only possible to have two physical solutions for P in the cases where 
a = 2, 3, 4. For these, the second solution becomes unphysical if c > 64. Thus only for the cases with 2 < a < 4 and 
Cmin < c < 64 do we get two possible static solutions to the Einstein equations for the same input parameters. For 
a = f , A = and the solution P_ = f is spurious; it does not solve the original unsquared eq. (^0)) . For a = 1, c is 
restricted to the interval [0, 64] for the valid solution P+ to exist, since P+ < for c > 64. For a = 5, Cmin = 64 and 
the solution P_ is always negative. For a > 5, although Cmin < 64, P_ is still negative and there is only one physical 
solution. We illustrate the three qualitatively different cases by graphing P versus c for a = 1, 3, 5 in figure [TJ 

We can also solve eq. (I^5t for the 3-brane tension in terms of P. To get a simpler-looking result, it is helpful to 
once again square both sides of the equation before solving. We find that 

^6^3 _ -, 5 /aT4fc|\ -5 ,^J^^^5+3 



2. =l-n^J (^^-1)-^-- (28) 

Recall that fc|f3 is the deficit angle of the conical singularity at the 3-brane, which cannot exceed 27r. Eq. I|28p is 
consistent with this requirement, so long as T4 > and P > f . But we should further demand that the brane tension 
be positive to avoid instabilities associated with negative tension branes. Eq. (|28p shows that this problem can always 
be avoided by choosing the free parameter T4, quantifying the amount of extra nontensional matter on the 4-brane, 
to be small enough. 
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FIG. 1: The solutions for P as a function of c = (k^iTi)^ for the cases of a = 1, 3, 5, where a determines equation of state of 
extra matter on the 4-brane. For 1 < a < 5 there are two solutions for P, as exemplified by the middle graph. 

III. MODIFIED FRIEDMANN EQUATION 

In this section we will use approximate and numerical methods to solve the jump conditions (j22l I23|) and thereby 
deduce the form of the Friedmann equation in the 6D model. Before solving them, it is important to understand 
which quantities should be considered as inputs and which are derived. Obviously we can freely specify all sources of 
stress energy, including Ag, T4, ra, T4 and a; these determine the Hubble parameter H. However our de Sitter brane 
solutions only exist for special values of P, the position of the 4-brane, which in the absence of the bulk scalar field is 
an unstabilized modulus, except in the case where a > 5 jTB]. Therefore eqs. ([22l [23|) should be seen as determining 
H and P given arbitrary sources of stress-energy. 

We will solve the jump conditions at first ignoring the bulk scalar field. The Friedmann equation is obtained by 
finding the dependence of H on the 3-brane tension T3 while holding other sources fixed. As shown above, for given 
values of the 4-brane stress energy, there exists a special value of the 3-brane tension (and possibly two such values) 
which leads to the static solution, H = 0. The analog of the Friedmann equation is the functional dependence of H 
on the difference St = T3 — which gives rise to expansion of the brane. Although this is not a realistic situation, 
since we are not allowing for radiation or matter on the 3-brane, one expects a pure tension source to nevertheless 
reveal the functional dependence of the Hubble expansion on the excess energy density of the brane. In any region of 
parameter space where the conventional Friedmann equation is recovered, this of course has to be the case, since the 
Friedmann equation depends only upon the total energy density and not upon the equation of state. To explicitly 
study brane matter with p ^ —p, it is necessary to regularize the 3-brane by giving it a finite thickness [l^, HEj. In the 
present work we avoid these complications by assuming that purely tensional energy density is sufficient for mapping 
out the functional dependence of H{p), but it is possible that this assumption could break down in the regions of 
interest, where significant deviations from the conventional Friedmann equation occur. This is an interesting question 
for future study, but beyond the scope of this paper. 

We will use two different methods to obtain the desired relation. First, in section UlI Al we will use a perturbative 
approach, treating H and 5t as small quantities; this gives some analytic insight into the deviations from the standard 
Friedmann equation. Second, in section llll C[ we solve the jump conditions numerically, allowing us to probe arbitrarily 
large values of the brane tension. 

A. Perturbative Friedmann equation 

A perturbative solution for H as a function of T3 can be found by expanding around the static solutions found in 
section IiTD] In that section we noted that for 1 < a < 5 there are two static solutions for a given set of stress-energies 
for the source brane. The procedure we describe applies equally to expanding around either of these two solutions 
in those cases. The strategy is to first solve eq. to obtain the position of the 4-brane in terms H, P = P{H). 
Substituting P into eq. (|23p. we then find a relation between H and T3 since the deficit parameter in eq. (j23p is 
given hy K = K{H,T3) from eq. (|2ip . The relation can be solved for H by perturbing in the excess 3-brane tension 
(5t3 = T3 — f3, or equivalently in the Hubble rate H. Thus we start by expressing P as a Taylor series in H using 
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eq. (|22p : then we will eliminate P from eq. ((23)) to obtain H as a Taylor series in St^. The Taylor expansion for P is 
given by 

- dP l(fP ^ 



where we have defined the dimensionless quantity 



n = H^f (30) 
The coefheient ^ in this expansion can be found by differentiating eq. ([22]) to get 

dP 4(a+ 1)P^ - 9a + 1 



dn fO(a-5)P5 + 5(3a + 5) 
To relate H to the 3-brane tension, we use (121]) and (l23l) to find that 



(31) 



'''^2^"lv 2pVV2G(P)^"-ir/2(i')y ^ ^ 

where we recall that p depends on H via eq. (|2T|) . Note that fa is defined to be the value of corresponding to 
H = and P = P. To find the Friedmann equation at leading order in St^, we take T3 = fa + St^ and substitute 
P = P + into the r.h.s. of This gives 

^ § ^-3 = -^W(P^ - 1) 1 - j + (33) 

which after solving for becomes 



i/^^^(P^-l)-^l-^r3j ^.^n^OiSr^) (34) 

Ignoring the terms 0(5t^), eq. ([M)) has the form of the standard Friedmann equation, = (87rG/3)p, since 5t^ is 
the excess energy density driving the cosmological expansion. However we need to check that the coefficient of 5t^ in 
eq. (|34| is indeed SttG/S. We can compute Newton's constant, or equivalently the 4D Planck mass = (SttG)"^, 
by dimensionally reducing the 6D Einstein-Hilbert action. Ignoring Kaluza-Klein excitations, the 6D and 4D Ricci 
scalars are related by Pg = Ri/a; thus 

Ml = {%ttG)-^ = j y/^a-^drd9 = j aK^/J^drdd 

= 2nM^Ki f dr = ^-^M^f^^ (l _ k^^) (35) 



2g^ 

Notice that in the above expression, the 4D Planck mass depends upon H through the integration limits g and P; 
thus Newton's constant becomes time-dependent when the rate of expansion of the universe is not constant. We will 
need to check whether this effect can be small enough to be consistent with experimental constraints on the time 
variation of G (see next subsection). But for the immediate purpose of comparing our obtained Friedmann equation 
with the standard one, we define a static Planck mass which corresponds to its value when H = 0: 

Ml = (87rG)-i = 2^MlKi ( dr ^ ^fM^{P^ - 1) f 1 - ^-^3 ) (36) 

Ji 15 \ 27r / 

Since M| — k^^, we see that indeed eq. ([M]) is consistent with the 4D Friedmann equation, — {8nG/3)ST3. This 
is in contrast to ref. [l8| . which mistakenly found deviations from the normal Friedmann equation at this order. The 
mistake made there was that corrections to the expansion rate due to changes in the 3-brane and 4-brane positions 
were argued to be negligible. However, we have just shown that it is essential to take account of the changes in g and 
P in order to obtain the correct result. 

We emphasize that this result holds equally in the vicinity of either of the two static solutions that exist in the 
cases of 1 < a < 5. This might at first seem contradictory since the slopes of versus St^ differ between the two 
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solutions. However once one fixes the value of the 4D Planck mass by appropriately choosing the values of the input 
stress-energies, either solution equally well reproduces the standard Friedmann equation at small values of H . 

Having established that the Friedmann equation is recovered at leading order in 5t^ , we now turn to the corrections 
at higher order. The detailed calculations are presented in the appendix. Using eq. ([5^ and solving for H^, we obtain 

= ^-^5r, (^1 + [^)'ji'Srl^ + 0(<5r|) (37) 

where J depends on P and a as shown in figure [21 At large P, J increases like P^ for generic values of a and like 
P^^ for a = 5 due to the behavior of the denominator in eq. (jA5|) . For intermediate values of P, and 1 < a < 5, J 
can become negative, as seen between the cusps of In | J| in figure [2] 




0.2 0.4 0.6 0.8 
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FIG. 2: log \J\ versus logP for several values of a. J is positive except for the regions between cusps. 



B. Time and spatial variation of Newton's constant 



From eq. ((35)) we see that the 4D Planck mass depends on H and t^, the 3-brane tension. We can therefore 
anticipate that it will depend on time in a situation where the pure tension T3 is replaced by matter or radiation 
which gets diluted by the expansion of the universe. To quantify this dependence, we will consider the variation of 
M4 (or equivalently Newton's constant G) with St^. 

For small St^ we can differentiate eq. psp with respect to ra, keeping in mind that P and g depend on 7i, which in 
turns depends on T3. At leading order, we find that 

d\nG _ k f 3 P^ {4.ia + l)P'^ -9a + l) 3 -3-3 ^ 



dST3 (p3_i)2^i_fc^3) 2 (a-5)p5 + 3a + 5 2 



P^P-" - 1) (38) 



at Sts = 0, where — fc|/(27r). Since the deficit angle should lie in the physical range [0, 27r), k^f^ must be in 
[0, 1) which means the energy scale determining the 3-brane tension must be approximately equal to or below the 6D 
Planck mass. In the limit of large warping, P 3> 1, this simplifies to 

c?(5t3 5 — a (1 — kf^) 

assuming that a 7^ 5. Let us compare this to experimental constraints on the time- variation of G, \d\nG/dt\ < 10"^^/y 
[2^. Since dp/dt = —2pot'^/t^ , taking po to be the present value '--^ (3 x 10"'^ eV)"* and t = to ^ 15 Gyr, we find that 
the constraint implies that 

dlnG ^ d^nGldt ^ 
doT-i dp/dt 

Therefore k^ ~ < 10'*'* GeV^'*, and we have only the extremely weak constraint that Mg > 10"^* GeV. 
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A stronger constraint can be obtained if we make the reasonable assumption that the 6D compactification volume 
and hence G would be locally perturbed by a high density of matter, similarly to the global perturbation of G due to 
the 3-brane tension. For example, G should not vary greatly from its known value in the vicinity of a neutron star, 
in order to maintain precision tests of general relativity from binary pulsars. From this we conclude that k should be 
less than the inverse density of nuclear matter, leading to the much more interesting constraint 

Mq > 1 GeV (41) 

For the special case a = 5, the above approximations are not valid, and the time variation of G is relatively large 
if the radial size of the extra dimension P is much greater than unity. We find that 

d\nG ^ 9-5 

~d6^^'5 (1 - fc2f3) ^ ' 

The constraint on Mq is more stringent in that case by a factor of P^. For the interesting value of P ^ 10^® which is 
required to solve the weak scale hierarchy problem, the constraint is fc^ < 10^"^^ GeV^^, which gives Mq > 10^ GeV, 
from the time variation of Newton's constant; the corresponding limit from the spatial variation is that Mg must be 
close to the Planck scale. However if the warping is weak, P ~ 1, we retain the previous constraint 

C. Numerical Results 

We can determine the relation between H and T3 numerically using a simple algorithm. Notice that from eq. (j22p , 
one can solve for the combination ^/x = \/ P"^ — P^^ + PP-(P-P in terms of P and other known quantities, since eq. 
([22|) is quadratic in ^Jx. Then z = (Hgl)'^ ~ x + 1/P^ — P^. Furthermore, we can solve for g as a function of z using 
eq. (PDl) . which can be written in the form g — {1 + z/g^)^^^^. Although this has no analytic solution, numerically it 
converges very quickly to the exact result, starting from the initial guess g = 1. At this point, we have determined 
H{P) since H = ^/z/^ig). The final step is to solve eq. (I23p for T3, which appears in the factor K through eq. (j21|l . 
Finally we have H{P) and T3(P), which allows us to plot the relation between H and t^. 

In figure [3] we plot H{t3) for the three illustrative cases a = 1,3,5 (recall that a determines the equation of state 
of the extra component of matter on the 4-brane, needed for obtaining a compact solution, with P < 00). When 
a = 1, the solution to eq. ([22|) allows for the position of the 4-brane P ^ while Hgi ~ P"^/^. On the other hand 
eq. ([28| shows that P ^ as the deficit angle approaches 27r. The result is that H diverges for finite T3, and this 
behavior is seen in figure[21 We graph the dimensionless quantity [HlY versus T-i/2n in units Mg = 1. In these units, 
t^/2tt is the deficit angle of the conical defect over 27r, which should therefore lie in the range [0, 1]. Notice that 
is linear in (5t3 in the vicinity of iJ = as expected from GR, and as we showed in the perturbative treatment. The 
departure from the standard 4D behavior becomes evident once the deficit angle starts to become large. We used the 
parameter values t = 79.925, T4 = 0.1, r4 = 0.01 in Mg = 1 units for the a = 1 graph. Recall from section Hi Dl that 
the parameter c = {k'^tT^Y' had to be within a rather narrow range of values to get the interesting solutions with two 
branches; this explains the special choice of 

For the cases 1 < a < 5, we showed in section Hi D I that there are two static solutions, and correspondingly, when 
H > there are also two solutions for y^. Hence we obtain two branches of the Friedmann equation, as shown in 
the second panel (a = 3 case) of figure [3l but very interestingly these two branches smoothly join with each other to 
make a single curve when plotted in the H-T3 plane. At small H, each of the two branches is a linear function and 
agrees with the standard 4D Friedmann equation, as we established in section Fill Al The two branches meet near a 
maximum allowed value of the expansion rate, i/max, whose value depends on a. This dependence is shown in figure 
m These examples use parameter values £ = 79.925, T4 = 0.1, T4 = 1. 

The third panel of figure[3]illustrates the behavior when a > 5. Similarly to the 1 < a < 5 cases there is a maximum 
value of H (close to, but not exactly at = 1/^), but unlike thoses cases, there is only a single branch which is 
connected to a static solution. The branch at small deficit angle asymptotes to nonzero H. The other parameter 
values used for this figure are the same as for a = 3. 

IV. FRIEDMANN EQUATION WITH GOLDBERGER-WISE SCALAR FIELD 

We have shown that the model under consideration can be consistent with cosmological constraints on the time 
and spatial variation of Newton's constant and the standard 4D Friedmann equation at low energies, even in the 
absence of a mechanism to stabilize the extra dimensions. However as noted in ref. [l^, there is an instability for 
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t^/2ji '^J'^'^ 



FIG. 3: {Htf' versus fractional deficit angle (proportional to 3-brane tension in 6D Planck units) for a — 1,3,5, and other 
parameters i = 79.925, Ta = 0.1, n = 0.01 (for q = 1), r4 = 1 (for a > 1). 




1.5 2 2.5 3 3.5 4 

a 



FIG. 4: Maximum Hubble rate (log(ffmax^)) as a function of a for the cases similar to the middle panel of fig. [3] 

a < 5 manifested by the radion having negative mass squared, which makes it important introduce a stabilization 
mechanism (and of course fifth- force constraints also forbid a massless radion). As in [l^, we adopt the method 
of Goldberger and Wise (GW), who introduced a bulk scalar field 0, [SCj, [SlJ in 5D models to stabilize the extra 
dimension. Our interest will be to see how this affects the cosmological expansion. 

The basic idea behind the GW mechanism is to demand a nontrivial background solution for a free bulk scalar field 
(j) by including potentials for (p on the branes which bound the bulk. These have the property of being minimized for 
nonvanishing valeus of 0. The competition between gradient and potential energy in the bulk causes the energy to 
be minimized at some intermediate value of of the radius of the extra dimension. In ref. [16] it was shown that the 
resulting radion mass is suppressed by an extra power (namely 0.75) of the warp factor relative to the 5D Randall- 
Sundrum model, due to the largeness of the azimuthal extra dimension at the Planck brane in the 6D model. The 
presence of this rather large extra dimension in fact makes the present model a kind of hybrid between the ADD js^l 
large extra dimension scenario and the warped one. 

We have given the equations of motion and boundary conditions for the scalar field already in sections III Al and 
IIIBI However for the numerical methods we will use in the present section, it is more convenient to transform to a 
different coordinate system, whose line element is: 

ds^ = M^{f)[-dt^ + e^"^dx^] + df^ + L^{r) d0^, (43) 

In these coordinates, the 3-brane always sits at the origin, and we can rescale the 3-brane coordinates so that 
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M(0) = 1. Defining u = M' /M, v = L' /L [1^, the Einstein and scalar field equations are 



2A 



Af2 



1 



2A 

qq: Auv + Qu' = ^ + kl{^cj,'^ -V 



(44) 



where we have defined A4 = iH^ . In the following analysis we chose a special form for the bulk scalar field potential 



V 



2fc6 



!,2 ,4 



32 



10 



(45) 



inspired by the method of [33|, in which a similar potential was used in order to find an exact static solution for 
the coupled scalar field and gravitational system. (Note that in 6D, (f) has dimensions of mass squared and h has 
dimensions of inverse mass.) In passing, we remark that when u — v, i.e., in the limit where a — hence = Sgg, 
an exact solution exists, given by 



1 



-W, </)' 



dW 



V 



~ 2 J 



W = 



(46) 



Here the bulk cosmological constant term Ag = — 10// fcg has been absorbed into V. In this singular case where a = 
the 4-brane has been pushed off to infinite radius and the extra dimensions are no longer compactified. Furthermore, 
the space does not smoothly close at r = 0, so the 3-brane would have to be replaced by a 4-brane with small radius 
in this solution. We will confine our attention to a > in the remainder. 
The boundary conditions at the 4-brane, located at f = P, are given by 



Vp + Ti 

2L(P)« 

0'(P) 



6 



u{P) 



v{P) 



v{P) ~ u{P) 
1 dVp 



(47) 
(48) 
(49) 



where Vp is the scalar potential on the 4-branc, which prevents 4>'{P) from vanishing, even with no such potential on 
the 3-brane. The simplest nontrivial choice, which we adopt, is 



Vp = -Xcj) 



(50) 



The boundary conditions at the 3-brane at f = are more subtle because of the vanishing of the r-r metric element 
there, L{Q) — 0, causing v = L' / L to diverge as 1/r. Nevertheless, the Einstein equations (|44p are well behaved as 
r — > because vu and r(f>' are finite, and v' + v'^ = 0. For the numerical computations, we deal with this complication 
by starting the integration slightly away from the origin, for example at rg = 10"^", taking 



u{ro) = ro(A4-V) 
M{ro) = 1 



(51) 

(52) 
(53) 

(54) 



while 0(0) (or (p{ro)) is to be determined. ^ 



^ Since only the ratio v = L' /L appears in the differential equation and not L or L' by itself, for numerical convenience one is free to 
rescale L and L' such that L(ro) = ro, and restore the factor (1 f^) afterwards. 
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The differential equations can be solved numerically using the shooting method. For a given value of H, we first 
make a guess for ^(0) and integrate the equations from f = until the value f = P, defined to be the point where the 
boundary condition ((47|l is satisfied. However the condition ((49|) will not generally be satisfied; one has to adjust the 
initial condition 0(0) to achieve this (or alternatively tune A given some value of (/'(O)). Using the Newton- Raphson 
method, the correct value of 0(0) can be found after several iterations. Once the solution is found, the corresponding 
3-brane tension can be calculated from eqs. ([5T| l48 |) . which imply 



1 



2tt 



2{v(P) ~ u{P))Lf 



(55) 



In this way we obtain the relation between H and T3 which we interpret as the Friedmann equation. 

Wc illustrate the results in figurc[5l for the cases of a — 1, 3, 4, 5, and other parameters given by 6 = 1, A = —4.19532 
(this was tuned so as to give 0(0) = 4 in the static case), T4 = 0.1, T4 = 4 x 10"^ and £ — 88. The shapes of the 
curves are qualitatively similar to their unstabilized counterparts in figure [31 but the range of possible values for the 
Hubble parameter is dramatically increased in the cases a > 1. Departure from the linear relation between and 
T3 becomes significant only at much larger values of H than for the unstabilized model, as we would expect. 




0.86 



0.1 



0.2 



— a = 5 
■ ■ ■ a = 4 

a = 3 



0.3 0.4 



0.5 



0.6 



FIG. 5; Friedmann equation with stabilization by GW scalar field for a = 1 (left) and a = 3, 4, 5 (right) 



To further demonstrate the effect of stabilization, we have computed the fractional deviation of from its normal 
value, as a function of the excess brane tension (5r3, for both the stabilized and unstabilized a = 5 cases. We define 
AiJ^ = iJ^ - ^^5t^, so the fractional deviation is 



AiJ2 ^ STrG^Ta 



3F2 



(56) 



which is plotted in figure [H One can see that the deviation is more than 100 times smaller in the stabilized system 
than the unstabilized one. For reference, we also show the deviation in for other values of a in the unstabilized 
model, in figure [71 We use different parameters to keep the 4- brane position P fixed for reasons that will be explained 
later. For a = 6 case the bulk is stable to some value of He while for other cases the bulk is unstable for all values of 
H which will be discussed in section IVl 

To elucidate the behavior of the GW field we plot 0(r) versus r for different values of the Hubble rate and the 
4-brane coupling A in figures [51 and [51 There are two qualitatively different behaviors depending on whether A is 
positive or negative; since the sign of 0'(r) at the 4-brane is governed by A, the two cases lead to either increasing 
or decreasing away from r — Q. Moreover the magnitude of changes very little over the bulk if |A| is large, e.g., 
A = —4.2, while it varies exponentially with r if A is small, e.g., A = 0.02. The figures also show that the radial size of 
the extra dimension is roughly inversely related to the magnitude of A. In both cases, the magnitude of is decreased 
when Hubble expansion is turned on. The dependence of 0(0) on H is shown in figures [TPl [TT] for the two choices of A. 
Although decreases with H in both cases, the system with large |A| is much stiffer and exhibits smaller fractional 
change in 0. 



13 



0.06 



0.04 



0.02 



without OW field 
with GW field x 100 




0.25 




FIG. 7: AH^/H^ versus ra for a = 1 — 6 in unstabilized 
FIG. 6: AH'^/H^ (defined in eq. ([56}) versus ra in the model, 
case with and without GW scalar field for a = 5 
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FIG. 8: (t>(r) for A = 0.02 and several values of H. 



FIG. 9: (b(r) for A = —4.2 and several values of H . 



V. FLUCTUATION ANALYSIS 

To gain further insight into the issue of stability of the extra dimensions, as it affects the cosmological expansion, we 
have undertaken a stability analysis, by perturbing around the background with arbitrary fluctuations of the metric 
and the GW scalar field. Such an analysis has already been done around the static solution in ref. [lB| ■ Here we repeat 
those calculations but now we perturb around solutions with nonzero Hubble parameter. Our goal is to quantify how 
the Hubble expansion affects the stability of the background solutions. We will consider both the unstabilized model, 
and the version which is stabilized by the Goldberger-Wise mechanism. 



A. Unstabilized model 



The unstable mode of our model is the radion, which is an admixture of fluctuations of the metric functions a, b, c 
in the ansatz 



ds^ = a{r, t) -df + e^"^5^jdx'dx' + &(r, t) d0^ + c(r, t) dr 



(57) 
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FIG. 10: 4>{Q) for different values of Hubble rate {Htf, FIG. 11: <^(0) for different values of Hubble rate {Hlf 
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(58) 



where the 4-brane stress-energy components in eq. ([T^ are now S^j^jy 
perturbation around the background solution is 



'Vogf^v and Seg = -Vgggg. The smah 



a{r,t) = ao{r){l - ai{r,t)) , Oq 

6(r,0 = &o(r-)(l-6i(r,t)) , bo = f{r)K\ f{r)^r^ 
c{r,t) = co(r)(l+ci(r,i)) , cq = /"^W 



(59) 
(60) 
(61) 



Using the ansatz di + 3Hdi — — m^ai and similarly for 6i and ci to expand the Einstein equations (tt) + (ii), (rr) 
and (tr) to first order we get 



2ai +bi-ci^O 
a, 
ao 



boj 



ao 



ao 



6;+3al + i(J^I-^^)6r + i(^+3^)c,=0 
2 \bo ao I 2 \bo aoj 



3a^ ^ 6^ 
2 ao bo 
1 fb'r 



ci 
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(3ai -F6i) -h6i?^— ai = 
ao 



(62) 
(63) 

(64) 



For the junction conditions at r = P, using Z2 symmetry across the brane for the (tt) and (99) components of Einstein 
equation, we find 



a' 6' 
a b 



= klVo^fc\ 



r=P 



a 



= klVgy/c\^^p. 



(65) 
(66) 
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The first order perturbation of the junction conditions is 



" 3a'i + b[ ■ 
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-^Veci-kl5Ve 



(67) 
(68) 
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where we used eqs. (|65l66p to simplify their appearance. Furthermore, since £4 oc ^h{P), we can expand Vb, Ve as 



^0 = ^4 + ^4-^4", 



ZDo 
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Ve=n + {1- ay^L^", SVe = -77^ «(1 - ") ^4 i^" 

ZOq 



By subtracting eq. from (|65p and using (|69l70p . we find 
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^4 



aT4kQy^ \ ao bo 



Then eqs. §7]) and §Sll lead to 
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Combining eqs. (|72|) and (I73|) . we obtain the junction condition at 4-brane 
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''1 Jr=P 



2 V^o ao 



(ci + abi) I 
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(69) 
(70) 

(71) 

(72) 
(73) 

(74) 



We note that there are not really two boundary conditions at the 4-brane; imposing the momentum constraint 
on the junction condition ([75)) leads to the same result as eq. (17^ . 

At the 3-brane, we impose the condition that the deficit angle is unchanged by perturbations around the static 
solution. If we consider a circle at radius r = p + e around the 3-brane, with circumference L and physical radius Z?, 
we therefore demand that L/ D be invariant in the limit that e ^ 0: 
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In this way, we obtain the boundary condition at 3-brane 

[foi + ci] I =0. 

L ^ -^J \r—Q 

It proves to be convenient to work with a certain linear combination of the perturbations, 

X = 3ai + 61 . 



(75) 



(76) 
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Using the components (i,t) -I- (i,?) and (i,r) of the Einstein equation, the variables ai, hi and ci can be expressed in 
terms of X, 
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(79) 
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where / denotes the derivative with respect to r. 

Now we are ready to analyze the radion mass spectrum. Combining cqs. 



we get 
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The boundary conditions ([71)) and (fTT)) become 
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(83) 



r=P 



Thus if we solve eq. (1801) under the boundary condition ([8ip and (|83p. we can obtain the radion mass spectrum. In 
the asymptotic region r — )■ oo, we have 
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(84) 



from which we obtain the simple relation between the radion mass and the solution of bulk equation of motion. 
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(85) 



r=P 



which is consistent with the result in ref. [T^ when H = Hj (l + ij2^2-j(i/5) _ Note that the radion mass is expressed 
in terms of P so to numerically compare results with different values of a P should be fixed as we did in section [TVl 
In ref. [Toj it was noticed that one can obtain an analytic approximation for the radion mass, using only the 
asymptotic behavior of the solutions for X, which is valid in the case of strong warping, P 3> 1. At large r, the 
differential equation for X simplifies to the form 



X" ^ -X' ^ 0{r-'')X = 



(86) 



which has solutions of the form X = ci + C2/r^. Therefore generically X'/X ^ 1/f'*, and we can ignore the terms 
proportional to X'/X in eq. (j85p . Applying this observation in the present case where we take into account the 
Hubble expansion leads to the result 



5(a-5) 



2(a + l)#2 j (l + 0(p-3)) 



(87) 



We can notice a correlation between the radion mass squared and the deviations from the normal Friedmann equation, 
parametrized by AH^/H^ in eq. ([5S)) . We find that \AH'^ / H"^] is larger the more negative is, as one would 
intuitively expect. In figure [7] we showed AH^/H^ versus 6T3 for a = 1 — 6. The smallest deviation, in magnitude, 
occurs for the largest value of a. 

The radion mass squared is always negative when a < 5, but for a > 5, it can be positive if the Hubble rate 
is below some critical value, defined by = If the warping is large, P ^ 1, and the relation 

H = H/{1 + if2^2^(i/5) ^ implies that H 9^ H; in that case the critical value of H is 



5(a-5) y^' 1 3/2 
2(a + l) ' 
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For the a = 5 case, where the anisotropy of the 4-brane stress tensor is provided by the Casimir effect, the leading 
contribution ai H — vanishes, and one must look to the subleading terms of order . Of course, the approximation 
([88| is only valid for large warping, P ^ 1. In the discussion that follows, we will use the exact value of He which 
comes from solving the radion mass eigenvalue problem numerically. 

In section IVTl we will make use of the interesting case of a > 5, for which the radion is stable even without the 
introduction of an external stabilization mechanism."^ For a > 5 we can pose the question of how much of the 
exotic Friedmann relation at high values of H is consistent with stability of the radion. To accurately address this 
issue, we have numerically solved the eigenvalue problem eq. (jSOp for the radion mass rather than relying upon the 
approximation ([87)) or ([88)) . The result is shown in figure [T2l which plots as a function of H'^, for a = 6, T4 = 0.1, 
£ = 88 and T4 = 1. The radion becomes destabilized when {Ht}^ ^ 0-5- In figure [T51 we show the Friedmann relation 
in the vicinity of the transition region between stable and unstable cases. 

It is interesting that the transition between stable and unstable radion occurs quite close to the joining of the two 
branches which make H double- valued. In figure [M] we plot the value of {Ht£)^ at the turning point as a function of 
{Hc£)^, with the curve parametrized by varying c = [k^TAt)^ . (We still use T4 — 1 here, but in fact the value of T4 
does not affect the result because it only appears on the right hand of eq. ([23)) and can be rescaled away.) Figure 
[TJ] shows that to within numerical accuracy. He and Ht coincide: thus the exotic branch of the Friedmann equation, 
which leads a double-valued Hubble rate, is unstable to decompactification of the radius of the extra dimensions. 

It is also interesting to know how strongly H deviates from the prediction of general relativity before the critical 
value He is reached. Let us define the standard Hubble rate as Hq = St^/SM^. We calculated H^/H^ as a function 
of c and show the result in figure 1151 There it is seen that the maximum deviation is no greater than approximately 
1.5 (although with extreme fine-tuning of c larger deviations might be possible). Such a moderate deviation might 
seem unimportant if it occurred early enough in the history of the universe. However we will show in section [VD that 
there can be interesting effects during inflation, due to the derivative of H with respect to becoming singular at 
the critical point. 




B. Fluctuations with Goldberger-Wise stabilization 



Carrying out the analogous steps as before, but now including the GW bulk scalar field, we derive the coupled 
equations of motion for the radion field X and the fluctuations of the GW fleld, 0i. Using the shorthand A'^ = a'^/a^. 



Unfortunately we are not aware of any physical models which give such an equation of state for the extra matter on the 4-brane. 
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FIG. 14: Hubble rate {Ht£f at the turning point (where piQ 15. H'^/Hi at the turning point for diflerent values 
the two branches of H come together) versus (HJ) , the where Hi = Sr^/SM^ is the standard Friedmann 

critical H above which the radion is destabilized. relation. 
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The boundary conditions for X are 
X' -AoX = 0\e 

/3(a + 3) ' 3(a-l) ' 3 6co 



(91) 



8 



-Bo- 



U2 i'2 



2 Bq ao -Bo 2i?Q 



'2t 



-(s; - a;) ((-7 + 3a) a; + (« - i)i?;) + - s;) - -m' - ^— (a; + i?;) - - B',)ky^ 

3^0^, 6coi/2 
2Bq ao -Do '^'P 



1 + 3a „' ' 7 — 3a 



OV'O 



OV^O 



-^0o' 
2Bn ° , 



2^6^! = 



X 



(92) 



In principle, the analysis is complicated by the mixing between the radion X and the scalar fluctuations 0i . However 
it was shown by numerical solution of the equations in ref. [iBl that this is a negligible effect, and so we can neglect 
the terms proportional to 0i and ip'i in eq. (|92p . Again, we can solve for the radion mass in the limit of large warping, 
r = P> 1, 
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The large-r behavior of the radion wave function X is altered by the GW field; numerically we obtain a solution which 
is shown in figure [TBI One observes that X'/X is negative, and that falls faster than 1/r^, but slower than 

l/r^ (figure fT7|) . The sign and magnitude of X' /X implies that the term proportional to (P'q in ([93|l is positive and 
that it is the dominant contribution which makes > in the absence of Hubble expansion. However, since the 
order of X'/X is always smaller than the second term in ([93)1 dominates when H becomes sufficiently large, and 
the radion becomes unstable. 
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FIG. 16: Radion wave function X and its derivative, with 
stabilization by GW mechanism. 
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FIG. 17: Log-log plot of X'/X versus r, whose slope 
indicates that X'/X falls faster than 1/r^ (but not faster 
than 1/r^) at large r. 



We can compare the results here to our earlier observations concerning the Friedmann relation in the stabilized 
system, figure [SJ There it is seen that the departure of from linear dependence starts at an a-independent value of 
H ~ 100/^. However in section lTVl we used large values of the GW field, whereas in the present section we treated it as 
a perturbation. To make a meaningful comparison, we have recomputed the Friedmann relation in the GW-stabilized 
system for smaller values of the GW field, by taking smaller values of the coupling A. Motivated by the results of the 
unstabilized case we investigate whether the radion mass also vanishes at the turning point of the Friedmann equation 
in this case. Taking A = 0.02 which corresponds to 0(0) = 1.53 x 10""'^^ ai H = 0, we solve the eigenvalue problem eq. 
()89II92|) and find the radion mass — —0.00037 and the deviation from the standard Friedmann equation AH^/H^ = 
61% at the turning point where the upper and lower branches join, i.e., (H£)'^ = 0.895. Although the radion mass 
is not exactly zero at the turning point, this could be due to our use of the unperturbed background solutions ao(r) 
and /o(7') from solutions without GW field. To make a precise comparison the exact background solutions with GW 
field would be needed, which would require further numerical investigation. 

In summary, the fluctuation analysis for the radion suggests that the exotic branch of the Friedmann equation is 
invalidated by the instability of the radion, even in the model with Goldberger-Wise stabilization. Moreover, we have 
not found examples where the magnitude of deviates from its GR prediction by more than 60% before the onset 
of the instability. 



VI. APPLICATION OF THE FRIEDMANN EQUATION TO INFLATION ON THE BRANE 

An important application of modified Friedmann equations which has been used in 5D models is to the study of 
inflation on the brane, beginning with ref. '2T]. There the simplest chaotic inflation is studied and the inflaton can be 
under 4D Planck mass Mp but still above the 5D scale M5. In that work, the high energy regime H ^ p was used to 
obtain novel results, like the enhancement of number of e-foldings and reduction in the scale of the inflaton fleld in 
chaotic inflation. Ref. [2^1 further showed that inflation can be sustained for steeper potentials than conventionally. 

We are not able to achieve such dramatic effects in the present model, because of our observation that the radion 
becomes unstable before the deviation from the standard Friedmann equation becomes very large. However, we will 
now show that interesting effects can nonetheless occur if inflation starts near the point of instability; namely, large 
running of the spectral index and breaking of the consistency relation between the tensor ratio and spectral index. 
To demonstrate these effects, we will use the modified Friedmann equation obtained in section Hill to study chaotic 
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inflation on the 3-brane with the potential 

V = ^toV^ (94) 

We assume that the inflaton field on the 3-brane is varying slowly in time and thus its energy density affects the 
cosmological expansion in the same way as the excess 3-brane tension in previous sections; it would be surprising if 
corrections to this assumption depended in a discontinuous way on the value of the equation of state w — p/ p. To 
achieve stability of the extra dimensions, we will consider the model with a = 6 and focus on its low-energy branch, 
which we have shown to be stable. 

We will denote the modified Friedmann equation by 

i/2 = HlT{p) (95) 

Hq — p/3AIp is the standard Friedmann equation and J-{p) quantifies the deviation of from the standard expres- 
sion. Notice that ^(0) = 1. For the braneworld model we are considering, with a = 6, we have found an analytic 
form for J- which provides an extremely good fit to the numerical results for the stable (low-energy) branch: 



HV) =\^l + fm-fm[l-\^^J j j (96) 

Here Vm is the maximum energy density as shown in figure 1131 The fit and the exact numerical result for is 
plotted in figure fTSl for a particular set of braneworld parameters, T4 = 0.00587, £ = 1500, and T4 = 4 x 10~^^, in 
units of Mq, which we shall also adopt to illustrate our findings in this section. The 4D Planck mass is given by 
Mp — 53.8 in Afg — 1 units for this background. We then find that is fit using the values — 0.33, p — 2.6 and 
Vm = 4 X lO^^Mp. We have chosen this value in order to make the tensor-to-scalar ratio close to the experimental 
upper limit, since one of our goals is to show that deviations in the consistency condition could be an observable 
signature of the braneworld model. 

On the other hand, the maximum value of is 1 -I- = 1.33 for these parameters; thus we do not get the "steep 
infiation" effect discussed by ref. [l^. However, this conclusion might change in the scenario where stabilization is 
achieved using the Goldberger-Wise mechanism. In the last section we showed that with a small GW field, T ~ 2.5 
at the turning point which is also the dividing point between the stable and unstable branches. With a greater GW 
field the maximum value of J-' will be even greater and thus it may be possible to achieve steep infiation. However in 
this section we will focus on the new effects which are unique to our 6D braneworld model. 
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FIG. 18: Exact numerical result (solid line) and analytic 
fit (eq. 1961 dashed line) to the Friedmann equation. 
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FIG. 19: Us — 1 versus A*" from eq. (|99p for different val- 
ues of Nm, defined in eq. (|109p . The case Nm = 00 is 
standard chaotic inflation. 
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A. Spectral index and consistency condition 

We now consider how the modification T affects infiation. For our purposes, the main effect is in the reiation 
between the power spectrum of tensor or scaiar perturbations and the inftationary potentiah 

P - -TP 

where Pg^ and Pt^o are the conventionaUy defined power spectra when expressed in terms of V ^ and e is the usual 
slow-roll parameter, 

1 ^ / 1/' \ ^ 2M2 

The factors of T introduce additional scale dependence in the power spectra beyond that contained in Psfi and Pt,o- 
This shows up in modifications to the spectral indices, 

Pt 16e , , 

n - ^ - ;^ (101) 

These can be derived using the horizon-crossing condition kjPl — a — and the relation between the number of 
e-foldings N and Lp, 

where the exact value of (T) is not crucial for our purposes, since only varies between 1 and 1.33; however we can 
estimate it as 

(J^) ^ — / TdV ^ 1.05 (103) 

''m Jo 

Henceforth we set {J-) to unity. These relations imply 

dlnfc = dAr+ -rflni/2 = ( _ + — ( i + \dLp = i^dV (104) 

In the slow-roll approximation, this is dominated by the first term, dlnfc = —(TVlM'^V'^^dV^ which allows us to 
carry out the differentiations leading to eqs. (|99m0ip . 

Before investigating the new effects due to T , we will impose observational constraints on the inflation model. For 
chaotic inflation, these are conveniently expressed in terms of the total number of e-foldings, 

where we used (/Sj = 2M2 by assuming inflation ends when e = 1. Therefore the initial value of the inflaton is given 
by ¥5-/Af2 = 4A^e- The COBE normahzation on P,. gives (H] 



2,„4 



eM4 4M6 



= A= 7.1 X 10"M(fco) (106) 
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where A(fco) = 0.7— 1.0 is amplitude at the pivot scale feg, which we assume corresponds to the field value Lpi. Eqs. 

p e 

Furtherfore eqs. (f98l llOli I105P show that the tensor-to-scalar ratio is given by 

n = (108) 

One consequence of the new dependence on T is that the spectral index depends differently on the number of e- 
foldings than in standard chaotic inflation, where — 1 = — 2/iVe. From eq. ((99|) . we obtain the modified dependence 
which is graphed in figure 1191 for different values of Vm , which we parameterize as 

V„r = ^mV™ = 2m^M^N„, (109) 

Here the parameter is the maximum number of e- foldings due to the limitation V < Vm, and we used the relation 
ip/Mp = AN between the field value and the number of e-foldings. Interestingly, there is an upper limit on 
depending on Nm which is significantly smaller than 1 — 2/Nm- 

A further novel feature is that the consistency condition between the tensor spectral index and the tensor-to-scalar 
ratio is different from the conventional prediction, 

Uf 1 dhiJ-\ , ^ 

Using the expression ((96)) for JF, we obtain sizeable deviations oirit/r from the standard prediction (—1/8) when the 
inflationary potential V starts out being moderately tuned to the maximum value Vm allowed by the brane model. 
The deviations as a function of V are shown in figures [20] and [21] 
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FIG. 20: 9,nt/rt versus V/Vm, showing strong deviations FIG. 21: 8nt/rt versus Iogi(,(l- V/V™), further illustrat- 
from the standard predicted value oi -1 as V ~* Vm- ing the deviation from the consistency condition when V 

is very close to Vm- 

The dependence of rit/rt on wave number in the CMB spectrum is the more observationally relevant issue. We can 
analytically elucidate this dependence using the fact that itself does not deviate greatly during inflation, even though 
its derivative is large. The approximation in eq. (|104p gives dhik/dV = — J- / {2m? M"^) , so that Infc is approximately 
linear in the potential itself. To demonstrate this, we plot the exact dependence of Vmdhik/dV, 

in figure [551 and the integrated result for In fc as a function of V -, for — 55 (for other values of the result 
looks very similar). The figure shows that Infc is to a very good approximation proportional to V . Therefore the 
wave-number dependence of rit jr can be inferred from figures [^Dj and 1211 
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FIG. 22: Vmdlnk/dV versus V/Vm, and the integrated 
result for Infc (plus an arbitrary constant, to put it on 
the same graph) versus V/Vm- Solid curve is Infc, dashed 
is the linear fit. 
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FIG. 23: logic as versus logio(f - N/N^) for Nm = 
30, 40, 50, 60 (top to bottom curves) in the m?if '^ chaotic 
inflation model. 



B. Running of spectral index 

Another possible spectral feature which can be enhanced by the braneworld modification is running of the spectral 
index, as — dus/dlnk. From eqs. (|104p and ((99|) we find 

1 / , M^V'V" dlnT ^dliiT yd^lnT ^fdliiTV\ , , 

- {^''^^ - ~ ^^Y^ + W - dh^ + W - Kd^ ) ) ^'''^ 

This reduces to the standard expression when J- = 1. But when the potential is close to its maximum value, as is 
enhanced by the large values of the derivatives of J-. This is potentially interesting because in standard inflation, as 
is of higher order in the slow roll parameters, and thus negligibly small, whereas observationally there is marginal 
evidence for running at the level of as ^ —0.1. 

However, in the chaotic inflation model, we find that as is still negligible unless V is extremely fine-tuned to be 
close to the maximum value Vm (equivalently, N must be very close to Nm)- As illustrated in figure [23l the required 
tuning is worse than 1 part in 10'' for N > 40, which includes the preferred values of N for the chaotic inflation model. 
Furthermore, eq. (jll2p predicts that as > due to the enhancement, whereas the marginal experimental indications 
are for negative values of the running. 



VII. CONCLUSION 



In this paper we have studied the warped codimension-two braneworld model which most closely resembles the 
original codimension-one warped scenario of Randall and Sundrum, with the goal of understanding how cosmological 
expansion (the Friedmann equation) is altered due to the extra dimensions. The bulk is approximately AdSg with 
a 3-brane at the bottom of the throat and a 4-brane at the top; in addition to tension, the 4-brane has an extra 
nontensional source whose pressure in the angular direction varies with the brane circumference ~ L4 as L^". This 
extra source is needed in order to find solutions with localized gravity when no stabilization mechanism is included, 
and it is known that the radion is stable when a > 5. It generally plays a less important role when stabilization 
is induced, for example, using the Goldberger-Wise (GW) mechanism. 

We determined the modified Friedmann equation both analytically, in a perturbative expansion in powers of the 
energy density on the brane, and numerically. The perturbative result, eq. p7p . has the curious feature that ~ 
p + 0{p^), i.e., the 0{p^) correction vanishes. For larger values of p, the numerical solution is necessary, as described 
in section UlI CI We checked that the two kinds of solutions agree with each other in the small-p region where both 
are valid. 

In our study of the modified Friedmann equation for this model, we have corrected a mistaken claim of ref. [18|, 
which did not recover the standard GR result at low energies. In fact we find standard cosmology at low energy 
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even for the unstabilized models. At high energies, we find exotic features, in which the Hubble rate H{p) turns 
over and joins with another branch of the function at some critical value of p and H, making H{p) double-valued — 
see figures |3] and [S]'' This behavior even persists in the model with GW stabilization, but we have found that the 
exotic branch is nevertheless unphysical, because the model becomes destabilized by the Hubble expansion at the 
turning point. We discovered this by doing a small fluctuation analysis around the de Sitter solutions in section fVl 
in which the eigenvalue problem for the radion mass was solved both using analytical approximations and exact 
numerical methods. This result is reassuring, since otherwise we would be left with the puzzle of why the cosmological 
background is not uniquely determined by the sources of stress energy in the Einstein equation. 

Nevertheless, we have found that interesting deviations from the standard Friedmann equation occur near the 
threshold for radion instability: dH^ /dp diverges at this point, even though H itself is finite. In fact for the specific 
example we considered, H only differs from its standard value by 30% at this point, but the divergence in dH^ /dp 
can have observable effects if infiation happened to start near the maximum value of the potential allowed by radion 
stability. We demonstrated such effects in the simplest model of chaotic infiation on the 3-brane. Figure [19] shows 
that the spectral index is lower than in standard chaotic infiation; for instance with Ne = 55 e-foldings, Us = 0.964 
for standard chaotic infiation, whereas it is 0.950 for the braneworld model with Nm — 65 as the maximum number 
of e-foldings (hence the maximum value of the potential is 65/55 times its value at horizon crossing). Thus without 
much fine tuning, the braneworld model can have a potentially measurable impact on the spectral index. 

We then explored the breaking of the consistency condition between the tensor-to-scalar ratio rt and the tensor 
spectral index n^, which may be observable if is close to its experimental upper limit. Figure shows that 8nj/rj 
differs by 50% from the standard prediction of —1 even when V is only within 20% of its maximum value, so again 
interesting signatures can occur without excessive fine tuning of the model parameters. 

It would be interesting to try to extend our results to a string theoretic realization of inflation in a warped throat, 
namely the Klebanov-Strassler background [s^ with stabilization by fluxes [3^. Although we have used GW rather 
than fiux stabilization in our model, ref. (37j has argued that an effective field theory description of flux stabilization 
might coincide with the GW picture. So far, efforts in string theoretic inflation with warped throats have focused 
on brane-antibrane 'ss'l or DBI inflation. However the effects we have discussed should also be applicable in 
string-based constructions with conventional inflation confined to a brane. The phenomenon of radion destabilization 
at sufficiently high Hubble rate should also occur in string compactifications, since the decompactified background 
is always a solution, and the radion, being conformally coupled, gets a mass correction of order in the early 
universe. The H^(jP' term in the radion effective action thus drives (j) away from its nontrivial minimum toward the 
decompactified (f) — Q vacuum for sufficiently large H , and it is likely that deviations from the usual Friedmann 
equation will be large near this point. Therefore it is possible that stringy effects due to the extra dimensions may 
be manifested during infiation even if the infiaton is confined to the standard model brane. 
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APPENDIX A: DETAILED DERIVATION OF MODIFIED FRIEDMAN EQUATION 

We need to determine the higher order terms in the expansion of P in powers of 7i. Again differentiating eq. 
we obtain 



d^P 
dW^ 



2P 



-160 -I- 32 - 288 a - 96 a^) P^^ + (-800 + 288 - 480 a - 32 a^) P^^ 
1224 a + 48 + 264 a^) P^^ + (-504 - 24 + 600 -I- 3000 a) P^° + (-633 - 549 - 159 a - 195) P^ 
-2160 a - 1056 + 144 a^) P^ + (-402 a + 230 + 306 + 378 a^) - 50 + 522 + 162 + 390 a 



25 (5-|-2P^a-10P^ + 3a) 



(Al) 



* The same phenomenon occurs in the similar model of ref. [23|. but it is not commented upon there. 
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Substituting P = P + + into the r.h.s. of ([32|) . as advocated in the text, we get the surprising result 



6T3 + 0{dTl) 



(A2) 



Namely, the second order correction vanishes. 

To find the leading correction, we need to go the higher order. Continuing the procedure, i.e., computing dPP/dH^ 
in (|^ . using eq. (P^ . we find an unwieldy result 



d^P 



(5f 2 - 6656 a* + 23552 



89600 a + 5f 20 - 64000) P^^ 



-58240 a-" 



156800 - 280000 + 8512 



56000a- 448a^)p2^ 

4\ a23 



-99840 + 106240 + 5760 + 630400 a + 64000 - 51200 a^) P 
54000 a + 43200 + 64800 - 32400 - 129600 a^) P^^ 
12880 - 655200 + 54880 + 490000 - 1680 + 1246000 a) P^° 
-194080 + 38160 - 94000 + 5040 - 570000 a - 495840 a^) P^^ 
-82350 a* - 47250 a - 180900 - 33750 + 76950 + 267300 a^) P^^ 
-47040 + 385280 - 280000 - 2016000 a - 336000 a^) P^^ 
35280 a"* + 170000 + 611360 + 19600 a - 10800 4 
- 145800 + 35100 a"* - 67500 a + 67500 + 167400 



485280 a^) P" 
- 56700 a^) P" 



1225000 a + 7560 + 35000 + 35280 + 1089200 - 98280 a^) P 
84700 a - 1620 - 111500 - 134460 - 136440 - 356040 a^) P^ 
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4050 - 33750 + 60750 a + 8100 + 12150 a'' 



51300 a^) P^ 



'18900 a* - 318500 a - 279720 + 17500 + 1 1340 - 558600 a^) P^ 
-9360 + 60048 + 13608 + 23000 + 56376 a^P^ - 12600 a) P^ 

+33453 + 5103 + 28875 a + 81270 + 85050 - 4375 



125(-10P^ + 2 aP^ + 5 + 3 a)' 



(A3) 



Substituting P = P+ j^n + 5^^^ + ^^W^ into the r.h.s. of 



we get 



where 



§5r3 = ^ ( 1 - %s ) (P3 - 1) - + 0(7i4)) 



(A4) 



J{a,n) 



(-1120 - 1056 a + 96 + 32 a^)P^^ + (5088 a + 600 + 144 + 3Ua^)P 
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-(3000 - 24 + 360 - 1800 a)P^^ 



-8 + 120 - 600 a + 1000)P^ 



-(-459 - 165 - 3087 - 2433 a)P^^ + (-108 - 4500 -f 900 - 900 a)P" 
-(-36 - 1500 + 300 - 300 a)pi° -f (185 + 1179 -f 783 - 99 a)P^ 
-(2250 a - 162 + 270 + 2250)P^ + (90 + 750 - 54 + 750 a)P^ 

(-405 - 81 - 375 - 675 a)P^ - 225 a - 125 - 27 - 135 a" 



25 ((2 a - 10)P^ + 5 + 3 a)^ (P^ - 1) 



(A5) 
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APPENDIX B: SYMBOLS USED IN THE PAPER 



TABLE I: Definitions of symbols used in this paper (A-M) 



Variable 


Description 


Deflnition or flrst appearance 


A 


amplitude of CMB fluctuations 


cq. (fTOSIl 




A 


COBE normalization parameter 


eq. ([TOBJ 






logarithmic derivative of ao , &o , co 


above en fl^^QH 




a, f 


metric components in r coordinate 


eq. dH) 




a, b, c 


metric components in stability analysis 


eq. (f57)l 




ao,bo, Co 


unperturbed solutions for a, b, c 


eqs. ([59||-((6T|) 




ai, fei, ci 


perturbations of a, b, c 


eqs. (|5gil-(f5T1) 




a 


eq. of state parameter for ta 


ea (IT31l 




as 


running of spectral index 


eq. (ITT2I1 




b 


parameter of bulk potential 


eq. (|45} 




c 


dimensionless 4-brane tension parameter 


c = (kltTA)^ 




St3 


deviation of the 3-brane tension from fs 


(5t3 — T3 — T3 






deviation of from standard value 


ea (1561 




^, V 


conventional slow-roll parameters 


e = M'^V''^ I2V^ n = 




(jij-v, gee 


induced metric on 4-brane 


ea (IT31l 




T 


modification factor for Friedmann eq. 


ea (l95ll 




G 


4D Newton's constant 


ea 




G 


4D Newton's constant in static solution 


eq. dSnj 




H 


proper time Hubble parameter 


above eq. (|20p 




Ho 


standard GR value of H 


in the last paragrph 


of |V Aland eq. ((95l 


Ht 


H at turning point in Friedmann eq. 


flg.rm 






critical H for radion stability 


fig.[T4l 




H 


rescaled Hubble parameter 






n 


dimensionless squared Hubble parameter 


H = H'^l^ 




J 


coefficient in modified Friedmann eq. 


ea (1371 




k 


wave number for CMB fluctuations 


eq. (HU 




ko 


pivot scale for CMB fluctuations 


eq. dTOnj 




kl 


6D gravitational constant 


kl = Mg"", eq. @ 




k^ 


reduced 6D gravitational constant 


k^ = kl/2-K 




K 


conical deflcit parameter 


eq. dS]) 




I 


curvature length scale 


f = -lO/Ae 




La, L{P) 


(4-brane circumference) /27r 


ea (IT3I1 




X 


coupling of GW fleld to 4-brane 


eq. dSnj 




A4 


4D cosmological constant 


A4 = 




As 


6D cosmological constant 


eq. ® 




m 


inflaton mass 


eq. dH 




rrir 


radion mass 


above eq. 




M,L 


metric components in f coordinate 


eq. dSJ 




Mi 


4D Planck mass 


eq. dSSll 




Ms 


5D mass scale 


eq. © 




Me 


6D mass scale 


Mg-* = kl 




Mp 


4D Planck mass 


Mp = M4 




Mi 


4D Planck mass in static solution 


eq. dMll 
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TABLE II: Definitions of symbols used 


in this paper (N-X) 


Variable Description 


~a 

Definition or first appearance 




number of e-foldings at horizon crossing 


eq. (tT05|) 


AT 


maximum number of e-foldings 


eq. (HDD) 


Us 


scalar spectral index 


eq. dM} 


rit 


tensor spectral index 


eq. (fTOOll 


P 


"i-Ulcllic JJUoiLiUil Hi UlilLb (Ji >L 


eq. (21 and below eq. ([19)1 


r> 

_r 


4-brane position for tlie static solution 


eq. (I21I and eq. ((23 




scalar power spectrum 


eq. (EZll 


p 


tensor power spectrum 


eq. (HZj 


<P 


Goldberger-Wise field 


in section HVl 


00 


unperturbed G\V field 


in section IV Bl 


<pl 


perturbation to GW field 


in section IV Bl 


"(^ 


inflaton 


eq. dMll 


V5i 


inflaton at horizon crossing 


eq. (IT02)) 


Vf 


inflaton at end of inflation 


eq. (fT05l) 


P 


energy density 


eq. (HSJ 


Q 


3-brane position in units of ^ 


eq. dU and below eq. ((T5)| 


r, f 


radial bulk coordinates 


eq. dS} and eq. (|43l) 


r t 


tensor to scalar ration 


eq. (IT08)) 




O- Ui drlit; LclibiUli 


eq. © 


fa 


3-brane tension for the static solution 


above eq. ([2i)) 


J4 


4-brane tension 


below eq. (|3]) 


T4, 


nontensional 4-brane energy density 


eq. (HSJ 


U, V 


logarithmic derivatives of , Ij 


u = M' /M, V = L' /L 


V 


inflaton potential 


eq. dMI 




maximum value of inflaton potential 


below eq. dM)) 


Vo 


4-brane energy density 


eq. dSnj 


Ve 


4-brane angular pressure 


eq. dZQj 


Vp 


bulk scalar potential on the 4-brane 


Vp = ~\<j> 


V 


bulk scalar potential 




w 


"superpotential" for bulk potential 


eq. dMJ 


X 


combination of metric perturbations 


eq. dZHJ 
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